N95- 32229 


TDA Progress Report 42-121 


May 15, 1995 


Towards Optimum Demodulation of Bandwidth- 
Limited and Low SNR Square-Wave 
Subcarrier Signals 

Y. Feria 

Communications Systems Research Section 
W. Hurd 

Radio Frequency and Microwave Subsystems Section 


The optimum phase detector is presented for tracking square-wave subcarriers 
that have been bandwidth limited to a finite number of harmonics. The phase detec- 
tor is optimum in the sense that the loop signal-to-noise ratio (SNR) is maximized 
and, hence, the rms phase tracking error is minimized. The optimum phase detector 
is easy to implement and achieves substantial improvement. Also presented are the 
optimum weights to combine the signals demodulated from each of the harmonics. 
The optimum weighting provides SNR improvement of 0.1 to 0.15 dB when the 
subcarrier loop SNR is low (15 dB) and the number of harmonics is high (8 to 16). 


I. Introduction 

This work was motivated by the need for near-optimum demodulation of the extremely weak signal 
received from the Galileo spacecraft. This demonstration is accomplished in the buffered telemetry 
demodulator (BTD). Since the BI D is a software demodulator, it is practical to tailor the processing 

more closed to the Galileo signal conditions than would be practical in other systems, such as the 
Block V Receiver. 

A limitation of the BTD is that the input signal has been recorded by the full spectrum recorder 
and contains only the first four harmonics of the originally transmitted square-wave subcarrier The 
subcarrier phase detector initially implemented in the BTD uses a windowing technique similar to that 
used in the Advanced Receiver II and the Block V Receiver [I] but modified for the four-harmonic case 
[ j. ere is a parameter, W SCi that is analogous to the fractional window width in a square- wave sub- 
carrier phase detector. As shown in Fig. 1, this phase detector results in a degradation (loss in symbol 
signal-to-noise ratio (SNR) due to harmonic truncation and phase tracking error), which does not mono- 
tomcally decrease as the number of harmonics is increased. 1 In fact, when the tracking error is large and 
when the harmonics are combined using the usual 1/n weighting for the nth harmonic, it is sometimes 
e ter to use only four harmonics than to use all harmonics. This suggests two things: First, it tells us that 
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the phase detector may not be using the harmonics optimally. Second, it indicates that the demodulated 
harmonics may not be optimally combined. 


The phase detector used in [3] is derived from a window used on a square-wave subcarrier loop. is 
phase detector may not be the optimum for a finite-harmonic subcarrier. As a previous work [2] indicates 
the higher harmonics get larger phase noise jitters. Therefore, the effective signal amplitude on the nth 
harmonic is no longer' 1/n but some number smaller than that. The optimum weights to combine the 
demodulated harmonics should account for the SNR losses due to the loop. 


II. Optimum Phase Detector 

Here we derive a phase detector (PD) that is optimum in the sense that the loop SNR is maximized. 
To show the derivation, let us first take a look at the current phase detector used in the BTEb The 
current phase detector is the product of the combined in-phase signals s/P d d k cos0 c (8/tt )^„ =0 tW n 
+ ll 2 jcos[(2n + 1 )d> 3C ] and the combined quadrature signals sfV d d k cos</> c (8/7r 2 ) £ n=0 w„(l/(2n 
+ l))sin[(2n 4- 1 )</>., c ] where the w n are the weights used to combine the quadrature signals and, in 

the current BTD, these weights are 


w n = 


sin[(2n + l)(7r/2 )W sc ] 
2n + 1 


The loop SNR using the current BTD is derived as 2 


Psc — 


aj3 2 
" yB sc 


Pd 

No 



1 

2E s /No 



where 

2 „ TsoUi personal communication, Communications Systems Research Section, Jet Propulsion Laboratory, Pasadena, 
California, October 1994. 
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where L is the total number of harmonics used in the phase detector, P d /N 0 is the data power-to-noise 
ratio, E s /Nq is the symbol SNR, and B 3C is the subcarrier loop bandwidth. 

Now in order to maximize the subcarrier loop SNR, Psc , let w k , k = 0, • • • , L - 1, be unknown and a 
be the same as before, and differentiate the loop SNR, p sc , with respect to Wk and set the expression to 
zero. We then have 


dpsc _ Z 2 ^ W k 1 P <i a 

dw k 7 2 B sc N 0 a + l/(2E s /N 0 ) 


= 0 


( 1 ) 


Since Pd/N 0 7^ 0, o ^ 0, and 7, D sc are finite, the above is zero if and only if 

7 - Pw k = 0 


That is, 


L - 1 L-l 

N W n ~ ^2 W n W k = 0 

72=0 72 = 0 


or 


L-l 

^ Wn(wn - Wk) = 0, for all k 

72=0 

which implies that 


= Wfc, for all n and k 


The conclusion is that the optimum weights to combine the quadrature signals in the phase detector 
are a constant for all (finite) harmonics. Note that, for infinite harmonics, the parameters /? and 7 do 
not converge; therefore, the above weights cannot be used for square waves. When the optimum weights 
are used in the phase detector, the loop SNR becomes 
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( 2 ) 


L Pd a 

Psc ~ B sc No a + l/(2E s /N 0 ) 


Using the optimum weights in the phase detector (called the optimum phase detector), we can improve 
the loop SNR by 9.5 dB over the current BTD with window size = 1, and by 1.1 dB over the current 
BTD with window size = 1/4 (see Fig. 2). The same figure also shows that, using the optimum phase 
detector, the loop SNR obtained by using only one harmonic is higher than that using the current BTD 
with the window size being either 1 or 1/2. Note that when we use only one harmonic in the optimum 
phase detector, we may still use all the available harmonics to demodulate the subcarrier. 



Fig. 2. Comparison in loop SNR using the optimum 
phase detector and the current BTD. 


Degradations due to a finite-harmonic subcarrier loop can be computed using the expressions given in 
[21 Degradations as a function of the number of harmonics are shown in Fig. 3 Clearly we can observe 
that, using the optimum phase detector, we obtained a lower degradation with more harmonics. This 

agrees with our intuition. 


With the increase of the loop SNR, that is, with the increase of the number of harmonics, the linear 
region shrinks. See the normalized S-curves shown in Fig. 4. As the number of harmonics approaches 
infinity, the linear region of the S-curve approaches zero. In other words, this optimum phase detector is 
only for a finite number of harmonics. 


III. Optimum Combining Weights in Demodulation 

The demodulated harmonics are currently combined with the weight 1/n for the nth harmonic. These 
weights are optimum if each of the harmonics of the subcarrier is demodulated with the same phase jitter. 
In our case, however, we know that if the first harmonic has a phase jitter with a variance of a , then 
the nth harmonic would have a variance of (na) 2 . The weight 1/n is no longer optimum. 
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Fig. 4. Normalized S-curves. 


To derive the optimum combining weights, we assume that the harmonics are combined using unknown 
weights b n . We then express the SNR in terms of the weights. Differentiating the SNR with respect to 
the weights and setting it to zero, we should obtain the optimum weights. 

The optimum weight to combine the demodulated (2 n + l)th harmonics is derived in the Appendix as 

, _ cos[(2n+ 1)0 ^J 1 

cos 4> sc 2 n + 1 

When cj) sc is assumed to have a Tikhonov distribution, 

cos(2n + 1)0 SC = /f^KV^cCos^] pn + l 1 

J nI 0 (p 9C /4) [ 2 0 “J d0 “ 

Assuming that we have 4, 8, and 16 harmonics, the degradations in symbol SNR versus the subcarrier 
loop SNR, using the optimum weights and the usual 1/n weights, are compared in Figs. 5 through 7. 

IV. Approximated Optimum Combining Weights in Demodulation 

Since the cosine function is “smooth” in the vicinity of zero, for small phase jitters, n(p sc , the expected 
value of cos (n(p sc ) can be approximated by 


£{cos(n0 sc )} % 1 - n 2 — 


The approximated optimum weights are 
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DEGRADATION, dB 




1 - (2 n + 1 ) 2 <r 2 1 

1 - a 2 j2 2n + 1 


( 5 ) 


Note that this approximation is valid only when 
for four harmonics, the symbol SNR degradation 


nifisc is small. Using the approximated optimum weights 
is only slightly more than that using the optimum weight 


as shown in Fig. 5. 
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V. Conclusion 

We presented an optimum way of tracking and demodulating a finite-harmonic subcarrier. We found 
an optimum phase detector in the sense that the loop SNR is maximized. The more harmonics used, the 
higher the loop SNR we obtain. However, the linear region of the phase error signal shrinks with the 
increase of the number of harmonics. Therefore, this optimum phase detector is only appropriate for a 
finite number of harmonics. Using the optimum phase detector, the loop SNR is about 9.5 dB higher 
than that of the current BTD using window size 1, and is about 1 dB higher than that of the current 
BTD with window size 1/4. 

For demodulation, we found the optimum combining weights that account for the losses due to the 
phase jitter. Compared to using the usual 1/n combining weights, the use of the optimum combining 
weights can improve the symbol SNR by 0.1 to 0.15 dB at a low loop SNR (15 dB) and a high number 
of harmonics (8 to 16). 
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Appendix 

Derivation of the Optimum Combining 
Weights in Demodulation 


After each of the harmonics of the subcarrier is demodulated, the signals from each harmonic demod- 
ulation need to be combined. Assume that the combining weight for the (2 n + l)th harmonic is 6„; the 
signal amplitude at the Ith symbol is 


s = 


2 ^ * 

\fPddi cos <p c — y b n 

^ 71=0 


1 

2 n 4 1 


cos[(2n + 1)(A SC ] 


(A-l) 


where P d is the data power, and <j) c and 4> sc are the phase offsets of the carrier and subcarrier, respectively. 
The noise variance is 


a 


2 




No r 

2 r^aym 


(A-2) 


Taking the ratio of the average signal power and the noise variance, we have the average symbol SNR of 
the combined signal: 


SNR 


£i£!i 

2 < 7 2 


£{(4/7T 2 )Prf COS 2 0c(En=O bn cos[(2n + l)0 3c ]/(2n + l)) 2 } 
~ Enid blN 0 R sym 


(A-3) 


Differentiating the symbol SNR with respect to 6 fc , k = 0, • • • , L - 1, we have 


d{SNR) P d cos 2 ^>c(4/7r 2 ) 


(Enld b ^oR S ym) 


„ V-J cos[(2n -I- l)0ac] cos[(2fc + l)<ft 3C ] 

2 i^ bn 2 E +1 2 fc + 1 


L— 1 




aym 


n=0 


71=1 
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cos[(2n 4- 1)0 5 
2n 4 1 
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(A-4) 


Simplifying the above equation, we have 


8 


cos[(2fc 4 1)0 sc] V - '' 1 i 2 \ s 7 

2k + 1 n " ~ o ” 

n=0 7i=0 


cos[(2n -1- 1)0 5 
2n 4 1 




(A-5) 
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Let fc = 0 and bo = 1; we have 


L - 1 L—l 

COS 4> sc 'Y^b 2 n — 2^ b n 
n=0 n=0 


cos[(2n + 1 )<t> 3c] 
2 n 4" 1 


= 0 


That is, 

L—i r 

y: b n COS (j) ac bn - 
71 = 0 

Finally, solving for 6 n , we have the optimum combining weights, 

k _ cos[(2n -f 1 )<fisc] 1 

n cos 4> sc 2n + 1 


cos[(2n -F 1)0 sc 
2n -F 1 


(A-6) 


(A-7) 


(A-8) 
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The Galileo low-gain antenna mission will be supported by a coding system that 
uses a (14,1/4) inner convolutional code concatenated with Reed-Solomon codes 
of four different redundancies. Decoding for this code is designed to proceed in 
four distinct stages of Viterbi decoding followed by Reed-Solomon decoding. In 
each successive stage, the Reed-Solomon decoder only tries to decode the highest 
redundancy codewords not yet decoded in previous stages, and the Viterbi decoder 
redecodes its data utilizing the known symbols from previously decoded Reed- 
Solomon codewords. 

A previous article [1] analyzed a two-stage decoding option that was not selected 
by Galileo. The present article analyzes the four-stage decoding scheme and derives 
the near-optimum set of redundancies selected for use by Galileo. The performance 
improvements reiative to one- and two-stage decoding systems are evaluated. 


I. Introduction 

This article is a follow-on to [1], which analyzed two enhanced decoding options planned for the Galileo 
low-gain antenna (LGA) mission: Reed-Solomon redecoding using erasure declarations and Viterbi re- 
decoding using Reed-Solomon corrected symbols. The analysis in [1] produced tables of gains achievable 
from enhanced decoding under an assumption of infinite interleaving for one, two, or four stages of Viterbi 
decoding, but no Reed-Solomon redecoding, and for one or two stages of Viterbi decoding, with or with- 
out Reed-Solomon redecoding, under the actual Galileo conditions of depth-8 interleaving. The present 
article looks at the case of four stages of Viterbi decoding and depth-8 interleaving. The four-stage coding 
system has been selected for implementation to support the Galileo LGA mission. 


II. Block Diagram of Coding Options 

A block diagram of the various coding options is shown in Fig. 1. A Reed-Solomon encoded data 
block is interleaved to depth 8 and then encoded by the (14,1/4) convolutional encoder. The Reed- 
Solomon codewords can have four different levels of redundancies, as depicted by the lightly shaded 
areas at the bottom of the code block in Fig. 1. The encoded data are modulated, passed over an ad- 
ditive white Gaussian noise (AWGN) channel, demodulated, and presented to a Viterbi decoder. After 



